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1. Introduction 

The Hamiltonian description of geophysical fluids has been extensively developed 
[1, 2, 3, 4] and is suited to analyse the conservation laws and to investigate approximation 
schemes which preserve the geometric invariants. These methods are also applicable to 
the Poisson-Vlasov equations for a collisionless plasma [5, 6]. Examples of Hamiltonian 
models of geophysical fluids are the canonical formulation of Rossby waves in a rotating 
sphere and the system of N vortices. 

Finite dimensional Hamiltonian systems are well understood in the integrable and 
uniformly hyperbolic limits, which correspond to an ordered quasiperiodic and a chaotic 
evolution [7]. The effect of small perturbations is amenable to analytical treatment in 
both situations, but the case of large perturbations and the transition from regular to 
chaotic regimes can only be explored numerically. In addition the effect of a weak noise 
and of round off errors in numerical computations deserves a special attention. The 
divergence of orbits and the memory loss rate of a given system are intimately related 
and reflect its dynamic behaviour. The Lyapunov exponents, and related indicators, 
specify the asymptotic separation rate of two initially close orbits [8], whereas the 
autocorrelation decay measures how fast the evolution loses its memory of a given initial 
condition. The asymptotic divergence of nearby orbits, the decay of correlations and the 
spectrum of Poincare recurrences are intimately related. In particular the correlation 
decay rate is related to the Lyapunov spectrum [9]. 

When a small deterministic or stochastic perturbation is introduced into the system, 
the perturbed orbit diverges from the unperturbed one (with the same initial condition) 
and we call forward error (FE) their distance at time t or the root mean square deviation 
at time t, if the perturbation is stochastic. The reversibility error (RE) is the distance, 
or the root mean square deviation, from the initial point, of its image after the evolution 
forward and backward for the same time interval t. The loss of memory with respect to 
the reference orbit is given, for instance, by the classical hdelity [10], hrst introduced in 
quantum systems [11]. 

For linear symplectic maps the asymptotic equivalence or RE and FE was proved 
in [12] and a relation with the hdelity decay was established. The hdelity behaviour 
for maps on the torus and the cylinder with random perturbations was hrst analysed 
in [13] and compared with the perturbation induced by round oh errors [14, 15]. The 
RE analysis was hrst introduced to investigate the global ehects of round oh errors in 
symplectic maps in [16]. 

In this note we examine the asymptotic behaviour of FE and RE for Hamiltonian 
hows with a small random perturbation. In the case of linear hows we recover, with 
simpler proofs, the results previously obtained for linear symplectic maps. If the phase 
space is a compact set such as the torus then the autocorrelation and hdelity can be 
computed and their decay rate is related to the growth of the forward and reversibility 
errors. 

Denoting with a{t) the FE and RE for a random perturbation, the asymptotic 
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decay of fidelity is given by F{t) For an integrable system defined on the 

cylinder we prove that the error growth follows a power law a{t) rs-/ et^ where /9 = 3/2 
if the system is anisochronons, /3 = 1/2 if it is isochronons. The exponent 3/2 is dne 
to the conditional mixing, known as hlamentation in the case of a rotating flnid. This 
kind of mixing is also responsible for the algebraic decay of the spectrnm of Poincare 
recnrrences [17]. These results easily extend to integrable systems of higher dimension. 
If the unperturbed system is not skew the result very likely still holds, as suggested by 
numerical computations, though the proof is no longer elementary. 

The proposed method could be applied to any hnite dimensional dynamical system 
with a small additive noise because it is based on the linear equation satished by the 
Gaussian stochastic process which dehnes the forward and reversibility errors. Inhnite 
dimensional extensions might be considered since the stochastic process satishes a linear 
partial differential equation. 

We have analysed numerically the behaviour of the anharmonic oscillator, the 
Henon Heiles model, and a system of point vortices [18, 19, 20, 21]. The results obtained 
for these models agree with the ones previously obtained for the standard map and 
the three body system [12, 15]. Except for the vortex model the Hamiltonian reads 
H = T{p) + V{q,t) — ex fit) where fit) is a white noise and a symplectic integrator 
was used, 4-th order accurate for the deterministic component, just as for the three 
body problem. For a discussion on numerical stochastic integration see [22]. In the 
integrable and chaotic regions the results proved for linear systems are conhrmed. In 
the transition regions the results might be interpreted following the model proposed for 
Poincare recurrences [17]. 

Random perturbations are introduced by the environment and by observations: the 
effect of observational noise in dynamical systems is actively investigated, see [23, 24] 
and references therein. When the system is integrable the effect of a single observation 
is small. However a sequence of observations may cause an exponential decrease of 
the correlations or the hdelity. The multiplicative noise changes the signal in a way 
dependent on the signal itself and may determine a faster decay of hdelity when the 
change of phase between two observations is large. 

2. Dynamical systems with additive noise 

In a previous paper we have considered prototype model maps dehned on the torus or the 
cylinder, whose invariant measure is the normalized Lebesgue measure. We consider hrst 
a ID dynamical system with a stochastic perturbation dehned by the Langevin equation 

dx ^_ 

— = ^ix) + efit) ( 1 ) 

where x G M and fit) is a white noise. We denote by ( ) average with respect to all the 
realizations of the stochastic process 


m) = 0 


mm) = dit -1') 


( 2 ) 
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We denote the deterministic solution with initial point xq with x = Sfixfi) and with 
f{x,t) the distribution function at time t corresponding to an initial value /(a;o,0) = 
f{xo). The distribution function / for the deterministic system satishes the continuity 
equation 

t) + ■^{^{x)f{x, t)) = 0 (3) 

The fundamental solution of the continuity equation is given by 


G{x, xq] t) =6{x- Sfixo)) = 


5{xo - S-t{x)) 


r]{xo,t) = 


dStjxo) 

dxo 


(4) 


viS-tix),t) 

since the solution x = Sfixo) is assumed to have a unique inverse Xq = S^fix). The 
solution for the distribution function f{x,t) corresponding to an initial distribution 
f{xo, 0) = f{xo) reads 


f{x,t) = 


r*+oo 


G{x,Xo-,t) f{xo)dxo = 


fiS-t{x)) 

r]{S-t{x),t) 


(5) 


When the noise is present we can still write a stochastic continuity equation given by 
(3) where <h is replaced by <h(a;) + ef{t). In this case after averaging over the noise the 
distribution function satishes the Fokker-Planck equation (see [10]). 




whose solution is given by 


fe{x,t) = 


^+oo 


Gfix,Xo;t) f{xo)dxo 


( 6 ) 


(7) 


We denote with Gfix,XQ]t) the fundamental solution of the Fokker-Planck equation, 
namely its solution with initial condition Gfix, xq; 0) = S(x — xq). In the limit e —)■ 0 we 
recover the fundamental solution of equation (3). We write the solution of the Langevin 
equation as x(t) = Se,t(xo) where 


Seyixo) = St(xo) + ex(t) = Sfixo) + eS(t) -h O(e^) 


( 8 ) 


The stochastic process E(t) is the limit of x(t) for e —)■ 0 and proves to be Gaussian 
since it satishes the linearised equation 




M{t) = ^'{St{xo)) 


Indeed S(t) can be written as a convolution of the white noise E{t) = f* K{t, s) f{s) ds. 
If <h(a;) is linear 

<h(a;) = u + fix (9) 

then x{t) and S(t) are equal. In this case the deterministic component is 

_ 1 


Sfixo) = e^^xo + (j){t) 


(pit) = UJ 


fi 


( 10 ) 
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Xit) = e 
Jo 


y(t-s) 


^{s)ds 


= e' 


e2/3i _ 1 


( 11 ) 


As a consequence the fundamental solution of the Fokker-Planck equation (6) reads 


exp 


Ge{x,t) = 


- {x)y 

2a‘^{t) 


{x{t)) = St{xo) = e^^XQ + 0(t) (12) 


A/27ra2(t) 

We notice that for e = 0 the fundamental solution of the Fokker-Planck reduces to 
5{x — Sfixfi)), namely the fundamental solution of the deterministic continuity equation 
(3). We remark also that the root mean square deviation for fi > has an exponential 
growth a{t) ~ ee^*. If /3 = 0 the motion becomes a translation and in this case 
a{t) = et^G_ 


3. Observables on the torus and correlations 

In order to have a chaotic behaviour an exponential divergence of the flow and the phase 
space compactness are required. This can be achieved by considering a dynamical system 
on the torus T. The torus can be seen as the interval [0,1], where the endpoints are 
identified. Any periodic function of period 1 deflned on M is a dynamic variable on T^. 
If f{x) is any function deflned on M we may construct a dynamic variable f'tix) in the 
following way 


+00 


/tW = X] 


(13) 


J = -CO 


provided that the decrease of f{x) at infinity is sufficiently fast to insure the convergence 
of the series. The function fj{x) is periodic of period 1 by construction and therefore 
can be expanded in a Fourier series 

+CXD 


Mx) = ^ fk' 

k=—oo 


,27r ikx 


(14) 


whose coefficients are given by 


fk = 


pi +00 „ 

-'0 -'0 


— 27r ikx 


f{x + j)dx 


j=-oo 


(15) 


+00 

E 

j=-oo 


rj+l r+oo 

—2TTik{x—j) f ( _ / „—2TTikx 


f{x) dx = 


f{x) dx 


If f{x) G L^(M) its Fourier transform f{k) for A; G M is defined by 

/■-l-OO 


f{k) = 


^—2nikx 


f{x) dx 


(16) 
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and consequently fk = f{k) for k E Z. 

If f{x,t) is the solution of the continuity equation on M the corresponding solution 
on T is given by 


+ 00 +CXD ^ + 00 

fj{x,t) = f{x + j,t)=^ / S{x + j - St{xo)) f{xo)dxo = 

■_ •_ J — oo 


]=-oo 


j=-oo 


(17) 


r*+oo 


Gf{x,Xo,t) f{xo) dxo 


By Gji^XjXojt) we denote the fundamental solution on whose Fourier expansion is 
given by 


+ 00 

GT{x,xo;t)= ^ 5{x + j - Sfixo)) 

j=-oo 


-\-QO 

^ ^2'n-ikx ^—27vik St{xo) 
k=—oo 


(18) 


By replacing (18) into (17) we obtain the Fourier expansion of f'^{x,f) according to 


h{,x,t) 


+ 00 

E 


k = — OQ 


fkit) 


^27r ikx 


fk(t) 





(19) 


We introduce also the following coefficients 

Mt)= r (20) 

Jo 

Notice that for t = 0 we have +(0) = fk- The new time dependent coefficients differ 
from the Fourier components unless the evolution is just a translation Sfix) = xo + ut 
or f{x) has support on [0,1] vanishing elsewhere. In this case 

fk{t) = fk{t) 


and the following equality holds 


hix) fri^x, t)dx= / h{xo) h{St{xo)) dxo 


( 21 ) 


3.1. Correlations 

We introduce two distinct dehnitions of the correlation on the torus according to 
C{t) = [ fr{x) fj{x, t)) dx - /2 = ^ f-kfkit) 

kj^O 

see equation (19) 


( 22 ) 


C(i) = / Mx,i) fj(S,(x,i)) dx,i - /|5 = W f-kft(t) 

kf^O 


( 23 ) 
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see equation (20). These definitions are equivalent when St is a translation or f{x) has 
support on [0,1]. In general we can write C{t) according to 

C{t) = Y. ^k',-kf-k Ak,,_k = f dxo (24) 

(k',k)^{0,0) 

When the vector field <h(a;) is linear, see equation (9), and the fiow is linear, see equation 
(10), then C(t) = C(t) for a sequence of times t = tm where tm = \ogm for m G Z. 
Indeed from equations (19), (10) and (16) we have 

/fc(t) = e-2"*'^^«/(A;e^*) (25) 


The coefficients fkifi), according to equations (23) and (24), are given by 


+ 0O 




(26) 


k' = — QO 


where Ay-k is given by 


e-i^k'-keP^) sin( 7 r(A;' - ke^^)/{Ti{k' - k e^*)) if ^ m 


Ak',-k = 


(27) 


5k',rak if = m 

As a consequence, for t = tm, from (26) and (27) we have 

A(«„) = A(«„) = 


(28) 


4. Noisy systems on the torus: correlations and fidelity 

For a noisy system the evolution of a given initial distribution is given by equation (7) 
and if the system is defined on the torus T the evolution of the corresponding distribution 
/t((c) defined by (13) is a periodic function in x which can be expanded in a Fourier 
series according to 

/ + CX) +00 

G^T{x,Xo]t) fo{xo)dxo= Y (29) 

■°° fc=-oo 

where GeT{x,xo-,t) denotes the fundamental solution of the Fokker-Planck equation on 
the torus and fj{x,t) the distribution function on the torus at time t 

+00 +O0 

G^T{,x,XQ]t) = Y Gfix + i,XQ,t) fi j(x,t)= Y + (30) 

j=—oo j=—oo 

Even in this case we introduce two distinct definitions of correlations as an obvious 
extension of the noiseless case 
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and 


Cfit) 



(/t( S^^fixo))) dxo - fo 


( 32 ) 


where 

U{t)= 

do 

The hdelity gives the correlation between a dynamic variable f{x) evaluated along the 
perturbed orbit at time t and f{x) evaluated along the unperturbed orbit at the same 
time. We have a unique dehnition of hdelity given by 


m 


hiSfixo)) (/t( Se,t{xo ))) dxo - /o = 

f.kh' r 

./ \ / //-» r\\ H 0 


(fc,fc')7^(0,0) 


(33) 


Jf-.l. Linear systems 

For a linear Langevin equation, namely for $(a;) = fix + ca, explicit expressions can be 
obtained. Indeed, to evaluate fe,k(t), we expand G^r{x,Xo;t) in Fourier series, taking 
(8) and (11) into account. The result is 


Ge,k{XO) t) 



dxe G^Y{,x,XQ]t) 



g —27r i kx 


exp 


f {x-St{xo)f \ 

V 2a^(t) ) 

y^2n a‘^{t) 


_ ^-2mk St(xQ) g-27r2o-2(t) 


The Fourier components of feT{x,t) are given by 

^+oo 


fe,kit) = e 




^-2nikSt{xo) ^ g-2^ a it) 


(34) 


(35) 


In order to evaluate fe,k{t) we recall that for a linear Langevin equation (1) where <F(a;) 
is given by (9) the solution Sey{xo) is given by (8) and (11) so that 


(e' 


-27rifcSe,t(a:o)^ _ ^-2Tr i k Stixo) / ^-2 tt i k exit)\ _ ^-2n i k St(xo) ^-2TT‘^k‘^<T' 


fie 


) = G 


As a consequence taking equation (16) into account we can write 

Afc(t) = Ut) 

The explicit expression for the correlation is 


cx) = J 2 f-tfM 


, —27r^/c^cr^(t) 


(36) 


(37) 


(38) 
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and 

(39) 

In this case the hdelity has the following expression 

Ffit)= (40) 

{fc,fc')7^(0,0) do 


5. Model systems 

We consider hrst three different prototype models 

i) Translations on the torns Th 

The held is <h(a;) = u and Sfixo) = Xq + ut and the Fourier coefficients are 

/tM = A(«) = e-"'““‘A (41) 


The correlations do not decay. 

When the noise is introduced the root mean square deviation with respect to the 
unperturbed trajectory is a = and the Fourier coefficients are given by 

fikit) = u{t) = (42) 

The decay of correlations and hdelity is exponential with decay time r = (27r^e^)“^. 
The diherence is that the correlations oscillate and decay, while the hdelity goes to zero 
without oscillations, due to the absence of the phase factor, 
ii) Anisochronous translations on the cylinder 

This is a skew system. Letting x = {x,y) G C = T x I be the phase space 
coordinates, # = (|/, 0) the vector held and S't(xo) = (xq + yot^Uo) fhe how, for an 
observable f{x) to which corresponds frix) on the torus we have 

hit) = Mt) ^ y (43) 


Even though fj does not depend on the torus label yo the coefficient fk(t) depend on 
it according to the previous equation. The correlation is computed by averaging on the 
interval I = [0, a] to which yQ belongs. Taking account of the normalization factor a~^ 
the correlation reads 


C(() = (?(() = 5^ lA 


^—2n ikyot 


dyo = \fk 


2 ikat 


k^O 


sm{n k at) 
n kat 


(44) 


as a consequence supposing that \fk\ < A/\k\ we have 

3^2 1 

<- 

vra t 

and the decay follows a power law. This slow decay is due to the conditional mixing, 
known in physics as the hlamentation phenomenon, which occurs because the frequency 
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on the tori has a continuous and monotonic variation. For an observable f{x, y) such that 
the average on x is constant, the result still holds provided that the Fourier coefficients 
fkiy) are and decay faster than l/y at inhnity. 

Any autonomous ID system in the region around a stable fixed point delimited 
by the separatrices is dehned, using action angle variables (0,j), by the vector field 
$ = (a;(j),0). If uj{j) is monotonic, introducing the new variables x = fi/iflTi) and 
y = u{j) we are back to the case considered above. 

When a stochastic perturbation is introduced the vector field becomes ^{x,y) = 
{y + ^^x(t),efy(t )) where ^x{y) and ^y(t) are two independent white noises. In this case 
the stochastic flow S,:^t{xo,yo) is given by 

X = xo + yot + + ewxit) y = yo + eWy(t) (46) 

where w{t) = f{s) ds denotes a Wiener noise and w^^\t) = w{s) ds. Recalling that 
{w‘l{t)) = t, {{wl^\t)fi) = f^/3, {wx{t)wy^\t)) = 0, the root mean square deviations for 
X and y are given by 

/ D\i/2 

ax = e(t+ — ] (^y = (47) 


The correlation on the cylinder for the observable frlx) is given by 


C(*) = 5^ lAt 


-^ikat sin(7rfcat) 
71 kat 


g—2-77^ k (t+F /3) 


(48) 


The hdelity has the same expression where the term within the square parenthesis is 
replaced by 1. If the frequency is the same on each torus, or if the y variable is not 
affected by noise, the decay of correlations or hdelity is much slower since the cubic term 
in the exponential is absent. 

iii) Expanding flow 

Given the vector held <F(a;) = fix, the flow is Sfixo) = c^^Xq and the Fourier 
coefficients are given by 

/ +00 

= fikeO') (49) 

■oo 


If f{x) G then f{k) G and consequently \f{k)\ decreases to 0 for \k\ —)■ oo. If 

\fk\ < ^ then ffike^fi decreases exponentially fast to zero with t and the following 

estimate holds 

OO 

\C{t) I < 2A^ ^ e-d* (50) 

k=l 

The Fourier coefficients entering the other definition of correlation are, for ^ N, 



+00 

/c' = —OO 




smTT{k' — ke^fi 
k' — ked^ 


( 51 ) 
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If = m e N then letting t^ = \ogm we have fk{tm) = fk{tm) = fmk where 
\fmk\ < A/\k\. If f{x) has a compact support on [0,1] the two dehnitions agree 

for any t G M. 

When the noise is present the root mean square deviation is given by 


a = 


e2/3t_ 



1/2 


(52) 


The correlations consequently are given by (38) and (39) , the hdelities by (40) where 
fkit) and fkifi) are given by equations (49) and (51) respectively. 

The explicit expressions of correlations for ^ N is 


a(i) = f-kfik 

kfkO 




and 


Cfit) = ^ f-kfk'e-''- - ’ -e 


{k,k')jk{0,0) 

The explicit expression for the hdelity reads 


i(k'-keP^) simrUk'- ke^^) 

7r(/c' — /ce^*) 


Ffit)= 

{k,k')^(0,0) 


^ sin7r((fc^-/c)e^* 


t: i{k’ — k)e 




For — m G N we have 


(53) 


(54) 


(55) 


Cfit^) = Cfit^) = Y, Ffitm) = Y IM' (56) 

kjkO kfko 


6. Reversibility error and related fidelity 

The reversibility error is the distance from the initial point of its forward evolution up 
to time t and backward for the same time interval. Any autonomous flow is reversible 
since the flow St has the group property S^t ° St = I. This property is lost when the 
system is stochastically perturbed. In this case the local errors e.^(±t) at times ±t are 
independent and we have S^_t ° S^^t I. 

We dehne a hdelity associated to the reversibility error according to 

Feuit) = f dxo frixo) {MSe-t o S^^fixo))) - /q (57) 

do 

The FE and the RE are related, and for linear hows we prove they are equal up to a 
factor ^/2. 

We consider hrst the translation perturbed with a white noise whose how is 
Se,tXo = Xq + ut + ew(t) so that S^^_t ° S^^tXo = Xq + e{w{t) + w{—t)). The reversibility 
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error is given by the root mean square deviation of this process (TR{t) = since 

{w(t)w{—t)) = 0. If the flow is expanding, namely = /3a; + ef(t) then we have 

SeyXo = xo +exnit) Xnit) = e~'^\{t) + x{-t) ( 58 ) 

As a consequence from 

XR{t) = [ - f{s-t))ds 

Jo 

we obtain the reversibility error 

1 _ p-2/3i 

(^R{t) = e{xUt)y^^ = ^ - ( 59 ) 

The reversibility error o'ji{t) is equal to the forward error a{t) (times ^/2) but for the 
system whose vector held is <h(a;) = —fix rather than <l>(a;) = fix. The forward how is 
expanding, the reverse how is contracting, but RE for the expanding how is equal to the 
FE for the contracting how. To recover the equivalence it is sufficient to consider the how 
on the torus dehned by the following noisy vector held $£ = {—fix+efx{t),fiy+efy{t)) 
where ^t(t) and ^y{f) are independent white noises. A simple computation shows that in 
this case (jR{t) = \/2a{t). The deterministic how in this case is area preserving. More 
generally linear Hamiltonian hows exhibit the same equivalence asymptotically. 


6.1. Quadratic Hamiltonians 

We consider the general case of a quadratic Hamiltonian HQX) = | x - Bx where x G 
for an arbitrary number of degrees of freedom d and B is a symmetric matrix. The 
equations of motion and the how are 


where the matrix is symplectic namely = J. When the system is 

stochastically perturbed the equations of motion and the stochastic how become 

dx H 

— = Ax+e^(t) A.^tXo = e^*xo + ex(t) x(t) = J f{.s) ds (61) 

The RE is computed starting from the reversed how 


S'.^.tS'.^tXo = xo+Xfi(t) 


and explicitly reads 


Xij(f) = e[e ^*x(f) + x(-f)] 


e * (^(s)-^(s-t)) ds 
(62) 


ds Tr 





^r(^) = ( \\XR{t)f > = 2e2 


■t 


(63) 
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( 64 ) 


It is not hard to prove that a{t) and aR{t)/\/2 are asymptotically equal taking into 
account the property of symplectic matrices. Indeed and have the same 
eigenvalues. 

The hdelity for a noisy reversed flow on the torus is dehned by 


Fnit) 



d^O /T2d (xo) fT2d {Se,-tSe,t^o) 



C^Xo /t 2d (xo 


(65) 


The reversed flow is given by equation (62) where the global error exnit) is dehned. 
Expanding fj 2 d{x.) in a Fourier series we hud 


Fnit) = ^ /-k/k' [ dxQi--- j 

(k,k') 7 ^( 0 , 0 ) do Jo 


dX0 2d ^g27rik'6Xfl(i)^ ^ 


/-k /k' e 


-47r2 k-E2(qk 


( 66 ) 


where 


(k,k')5^{0,0) 


—As „—A^s 


S"^(t) = e ds e ^ e 

Jo 


(67) 


We recall that the square of the reversibility error is expressed by 


4(t)=6^Tr(S2(t)) (68) 

If all the eigenvalues are complex of unit modulus then is conjugated to an orthogonal 
matrix. If is orthogonal o'^(t) = 2de^t supposing the phase space is If A has 
real eigenvalues ±Aj ordered in an increasing sequence —A, —A 2 ,..., A 2 , A then is 
conjugated to the diagonal matrix = diag(e^*, ..., e“^*). Supposing for 

simplicity that A = A is diagonal (in the general case only trivial constant factors 
appear) we obtain that asymptotically 

(T(t) ~ (TR(t) ~ V2ee^^ k ■ S^(t)k =(69) 


Also in this more general case the asymptotic decay rate of hdelity is related to the 
asymptotic growth of the error. 


7. Models 

We give some examples of application of the previous results. To this end we consider 
as integrable system the anharmonic oscillator and as non integrable system the Henon 
Heiles Hamiltonian. The Hamilton’s equations are integrated by using symplectic 4- 
th order integrators, based on the splitting method. Finally we consider the system 
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of vortices for which explicit symplectic integrators are not available: we use then a 
Runge-Kutta integrator. The restricted 3 body problem has been discussed elsewhere 
[12] just as native discrete systems expressed by symplectic maps [14], 

Denoting the time step by At = T/ug, where T is a characteristic time of the 
system, and using fourth order methods, we insure that the local integration error is 
close to the round off choosing Ug = 10'^. In addition symplectic integrators insure that 
the error on the invariants has a null average growth. We consider also the Lyapunov 
error dehned as the distance at time t from the reference orbit of an orbit whose initial 
point has an error of amplitude e. We have checked that for integrable systems the 
forward error due to a stochastic perturbation of amplitude e grows as a{t) ef3/2 just 
as the reversibility error up to a factor a/2. However when the system is isochronous the 
growth is a{t) ~ We have checked that decay of correlations and hdelity follows 

an exponential law according to C(t) F{t) e integrable systems 

in the regions of regular motion, the error growth follows a power law with exponent 
3/2, whereas in the chaotic regions the growth is exponential a{t) ~ ee^*, where A is the 
maximum Lyapunov exponent. 

7.1. The anharmonic oscillator 
The Hamiltonian is given by 

Hg = — epf{t) = uj + r] cos'^6' + e(2j)^/^ sin6' f{t) (70) 

The transformation to new action angle variables (0, J) allows to render H integrable 
up to a remainder of order and neglecting this remainder as well as terms of order 
Ae the Hamiltonian reads 

H = uJ + + e{2jy^‘^ sinQ f{t) (71) 

In these new coordinates the noise become multiplicative. A much simpler model 
to solve analytically is the one in which the noise is additive 

H = HoiJ) -Qeej 0(t) + J e ee feit) HoiJ)=uJ + vl^ (72) 

where e = (ee,ej) is a unit vector and ^e, ^j(t) are independent white noises. 

In this case the equations of motion can be solved explicitly by retaining only the 
hrst order in e and for initial conditions J(0) = Jq, 0(0) = 0o read 

Q = Qo + il{Jo)t + €^l\Jo)ejWj\t) +€eewe{t) J{t) = Jq + €ejwj{t) (73) 

where D = dH^/dJ and D' = dQ/dJ . As long as the distance on the cylinder agrees 
with the Euclidean distance we can write for the forward error 
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1 

Figure 1. Plot of the exponent a for the fit a{t) = C en°‘ to the forward error of the 
Enharmonic oscillator with a stochastic perturbation of amplitude e as a function of 
the nonlinearity strength y. The noise is additive in Cartesian coordinates, according 
to equation (70). The red line is drawn just to guide the eye. When 77 <C 1 the system 
is almost isochronous and a = 1/2, when 77 3> 1 the system is strongly anisochronous 
and a = 3/2, but in the intermediate region the exponent varies continuously between 
these two values. 

The figure shows the result of a fit to the exponent a for power law growth of the forward 
error a{t) = Cef" as a function of the anisochronicity strength rj. The function 0 ( 77 ) 
varies between 1/2 and 3/2 with a rather sharp transition in agreement with theoretical 
estimate (74). 

1.2. The Henon Heiles Hamiltonian 

This model describes the motion of an elliptical galaxy and the Hamiltonian reads 

^ {pI +pI) + \+ y‘^)+x^y± ^y^ (75) 

The Hamiltonian i7_ is not integrable, whereas the model i7+ is integrable. For the 
first one we have considered the reversibility error due to the round off. In Figure 2 we 
compare the colour plot for the reversibility error with the phase portrait. Both refer 
to the Poincare section y = 0,y > 0 

7. 3. The N vortex model 

For inviscid plane fluid the dynamics of point vortices provides a complete description 
of the dynamic evolution. The Hamiltonian for N point vortices located at Zk = Xk + iyk 
with a strength Fj, is 


TfW = -fi ^ r,rsiog|2, - zt 


( 76 ) 
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Figure 2. Henon-Heiles non integrable Hamiltonian i7_ given by equation (75). Left 
side: colour plot of the reversibility error due to round off for H- = 1/8 and initial 
conditions on a grid of points in the (x,px), with y = 0,py >0 and t = 100 dynamic 
periods T = 27r. Right figure: phase space plot for the projection into the {x^px) phase 
plane of the orbits in the Poincare section y = 0,py > 0 
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Figure 3. Level lines and isocolor for the three vortices Hamiltonian in the symmetric 
case A = 1/2 


The coordinates Xk and yt are canonically conjugated. 

There are two independent hrst integrals in involution. As a consequence the system 
of 3 vortices is integrable, while a system of 4 is not. If the vorticity of one of them is very 
low, a reduced Hamiltonian can be written which looks similar to the restricted planar 
three body problem. Introducing a canonical change of coordinates the Hamiltonian 
reads 

H = - 


+ (1 - A) log[(a; + A)^ + log[(2^ + A - 1)^ + (77) 
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We remark that the splitting cannot be used to construct an explicit symplectic 
integrator. Only symplectic implicit schemes are available. However one can use a 4-th 
order Runge-Kutta integrator choosing Ug = 10"^ such that the Hamiltonian is conserved 
within an error very small with respect to the stochastic perturbation. In Figure 3 we 
show level curves of the Hamiltonian (77). The behaviour of the forward and reversibility 
error follows the power law with exponent 3/2 in agreement with the theory. 


8. The observational noise 


To complete the present analysis on the effect of noise on a dynamical system we consider 
now the effect of observations. We suppose they are made regularly with a time interval 
T. We assume our system to be regular and to be given by the translation on the torus 
T and that each observation, having a small duration r T, behaves as a white noise 
of amplitude e. The correlation and the hdelity for an observable fT{x) during the hrst 
measurement T < t < T + r decay. After this time interval the correlation oscillates, 
whereas the hdelity remains constant, namely for t > T + r is 

Cfit) = J2 (78) 

kJ^0 kf^0 


If we repeat measurements of the same duration r and same time separation T, after n 
of them, namely for t > n(T + r), the decrease factor becomes signihcant if 

ne^r > 1. The threshold to observe the exponential decrease is 


T 

r 


(79) 


8.1. Map formulation 

If the disturbance introduced by the measurement is proportional to the signal itself, 
then we must treat it as a multiplicative noise. In the continuous time case analytical 
estimates are difficult to obtain also for translations on the torus. To overcome the 
difficulty we consider the case of a discrete time system. When the noise is additive 
the results are similar to the continuous time case. We consider a map M which gives 
the evolution on the same time interval T which separates two subsequent observations. 
Letting Xn+i = M{xn) we have 

Xn = Xn-i +uT + e fn mod 1 (80) 
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where fn are random variables uniformly distributed in [—1,1]. The correlation after 
n = 1 iteration (observation) is 

Ce(l) = [ dxo hixo) {h{xo + uT + e = 


fkf-k’ e 

A:,fcV(0,0) 


271 ik ljT 




'-1 


(81) 


El/' 

k^ 


k\^ S{eT^/^ k) 


S{x) = 


sin(27rx) 

271X 


After n iterations (observations) the correlation becomes 

Cfin) = |/fc|2 S^{eT^/^k) (82) 

k^ 

The same result holds for the fidelity where the phase factor is absent. We recall that 
\S{x)\ < e“^^'°g(27r) i^i ^ ^ |5'(a;)| < (27r|a;|)“^ for |x| > 1, see reference [13], and 

follow the same procedure outlined there. Assuming -C 1 and that the Fourier 

coefficients are bounded by \fk\ < we have 

|C'(n)| < (83) 


8.2. Multiplicative noise 

If the random perturbation introduced by an observation is proportional to the signal 
effect and we make a single observation after N iterations the map (82) is replaced by 

ain = (x„_i+a;T)(l + modi (84) 


After the first iteration the correlation reads 


cfii)= Y 

k,k'^{0,0) 


277 i ku)T 


^^^2-KikujT g27ri(A:(l+eTl/2^)-fc')xo _ 




—k' ^77 i kcjT 


E „ 

( 0 , 0 ) 

2niker^^'^f (^T + 2 ^^ 


dfS 


A:,fcV(0,0) 

■ exp 




k — k' + \ 




(85) 

If f^ix) is a trigonometric polynomial, namely fk = 0 for \k\ > K and 1 then 

S{{k -k' - A;erV2^)/2) ~ 4,fc' so that 


c.(i) = E1/ 

k^O 


2^2.ikru.Tgf j^^^l/2 






u,T+l 


( 86 ) 


Supposing that uT 3> 1 a single observation produces a significant decrease in the 
correlation roughly proportional to (r^/^ ecnT)”^. On the contrary if 1 the 
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decrease is negligible and we need to make several observations. After n observations 
we have 


Cfin) = (^a;T+ (87) 

Supposing that 1 and uT 3> 1 we have the following estimates 

{ exp (- log(27r) nre'^iuTfi) if caT < 1 

( 88 ) 

{2TTT^^^eujT)-^ if r^/^ewT > 1 

To conclude, if the stochastic perturbation, introduced by an observation, is 
multiplicative rather than additive, a phase factor uT is introduced in the correlation 
decay. Since the perturbation amplitude and duration are assumed to be small 
(r^/^e 1), if the phase factor is large {uT 3> 1), the decay is much faster. If 

~ 1, a few observations determine a signihcant decay of correlations in the 
multiplicative case. 

9. Conclusions 

We have examined the asymptotic behaviour of the forward and reversibility errors (FE, 
RE) introduced by stochastic perturbations in a dynamical system dehned on a compact 
phase space and we have compared these errors with the decay rates of correlation and 
hdelity. We have first considered the translations on the torus T and the cylinder, 
showing that both errors follow the same power law with an exponent equal to 1/2 
and 3/2 respectively. For expanding flows on T the FE grows exponentially, but the 
equivalence of FE and RE occurs only if we have a linear flow on the torus which 
preserves the area. This is a special case of linear Hamiltonian flows for which the 
asymptotic equivalence of FE and RE is proved. We have analysed the relation between 
the memory loss and the error growth. For linear Hamiltonian systems we have shown 
that the correlation and hdelity decay as where a(t) is the FE and that a 

similar relation holds between the hdelity and the RE. The RE method can be used 
to investigate the transition regions from regular to chaotic motions and the ehect of 
round oh errors in numerical simulations. The proposed framework is general and an 
extension to classical helds, by using the theory of linear stochastic partial diherential 
equations to dehne the FE and RE, could be considered. Finally we have examined the 
ehect of the observational noise showing that for an integrable system only a sequence 
of observations can lead to the decay of the correlation and that the multiplicative noise 
can increase the decay rate if the phase advance between two observations is large. 
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Appendix A. Correlation on the torus and Perron Frobenius operator 

Given a map M{x) continuous on the torus there is a unique dehnition of correlation. 
Indeed the evolution of a given initial distribution fjlx) is governed by the Perron 
Frobenius equation 


fT{x,n) = {P^fT){x) = 


^ S{x- M^{xo))h{xo)dxo = ^ 


Xi 


(Al) 


where P is the Perron-Frobenius operator. Supposing M{xq) has q preimages, with 
Xj = {M~'^)j{x) for 1 < j < we have indicated all the preimages of M'^{xq). For any 
A C the invariance of the probability measure /i is expressed by /r(M“^(A)) = fr{A). 
For example the Bernoulli map M(x) = qx mod 1 where g G N has g inverses and 
M~^ = UkMfi^{A) where yi{M^^{A)) = q~^yi{A). Since the density of the absolutely 
continuous invariant measure is 1, we dehne the correlation according to 

C{n) = [ h{xo)fT{M'^{xo)) dxo - /q = 


dxo hixo) / 5{x- M'^{xo)) h{x) dx - = 


(A2) 


Mx)fT{x,n)dx - /o 


where the result has been obtained by interchanging the integration order. 

As a specihc example we consider the linear map M{x) = qx mod 1 continuous on 
T for integer g. In this case the preimages Xj of are 




_ ^ x + j 

n—1 qu qfi 


(.43) 


where jj = 0,1,..., g — 1 and j = ji + gj 2 + • • • + O'” ^ jn-i + (f" ^jn has a range from 0 
to g*^ — 1. As a consequence, from the Perron-Frobenius equation we have 


q"-l 


Mx,n) = g ”/T(g ""(x + j) 

j=0 


(A4) 



Errors, correlations, fidelity for noisy Hamiltonian flows. 
and the following equality can be proved 


21 


[ Mx,n) frix)dx = [ — fj (f^{x) = 

Jo r. Jo J \ J / 

(dl5) 

/-O+l)/?" /•! ^ ^ 

= Y dxo fj{xo) h {q''Xo - j) = / dxofrixo) frig'^xo) 

j=0 J3/<1" do 

where we have made the coordinates change x + j = q'^Xq. Letting be the Fourier 
coefficients of fj^x) the order k Fourier coefficients of fj{x,n) are fk{n) = fq^k- The 
equality (A4) is a consequence of continuity of the map M on T. For a linear flow 
Sfixo) = xq we have seen that a similar equality, given by equation (22), holds 
only for the sequence of values of tm such that = m is an integer, since the map 
^tmi^o) = mxomodl is continuous on T. 
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